A more accurate treatment of the problem of 
drawing the shortest hne on a surface* 



Leonhard Euler 



1. For a surface on which a shortest hne is to be drawn, given this differential 
equation between the three orthogonal coordinates x,y,z: dz — fdx + gdy, 
where / and g are functions of both x and y, let it be 

df = adx + l3dy and dg = f3dx + "fdy. 

Having done this, since an element of any line drawn on this surface would be 



^/dx^l^dy^l^, 

with the previous value put in place of dz an element of this curve will be 

= ■\/ dx'^ + dy'^ + {fdx + gdyY; 



whence if we put dy = pdx, this element will be dxy^l + pp + (f + gp)^ . 
2. This integral formula, which is to be brought to a minimum, will be 



dx^l+pp+ {f + gpf, 

which I indicated in general by / Zdx in my Treatise Methodus inveniendi lineas 
curvas Maximi Minimive proprietate gaudentes, so that in this case it would be 



Z=^l+pp+{f + gpf. 

Then, having put dZ = Mdx + Ndy + Pdp, I have shown that the nature 
of a Minimum or a Maximum is such that it is expressed by this equation^ 
Ndx = dP, which clearly leads to differentials of the second order. 



'Presented to the St. Petersburg Academy on January 25, 1779. Originally published as 
Accuratior evolutio problematis de linea brevissima in superficie quacunque ducenda, Nova 
acta academiae scientiarum Petropolitanae 15 (1806), 44-54. E727 in the Enestrom index. 
Translated from the Latin by Jordan Bell, Department of Mathematics, University of Toronto, 
Toronto, Canada. Email: Jordan. bell@utoronto.ca 

^Translator: This is namely the Euler-Lagrange equation in the calculus of variations. 
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3. Therefore since Z"^ = 1 + pp + {f + gp)"^, let us differentiate this formula 
and let the elements be separated into the three types, namely dx, dy, dp, and 
then be represented in this way: 

ZdZ = dx{a + f3p){f + gp) + dy(/3 + 7p)(/ + gp) + dp{p + g{f + gp)). 

Then since I have put in general dZ = Mdx + Ndy + Pdp, we will have in this 
case: 

^ ^ {a + I3p){f + gp) 
Z 

+ ip)if + gp) 



N = 

P = 



Z 

P + 9if + 9P) 



Z 

Hence (since (5dx + ^pdx = dg) it becomes 

Ndx = ^iil+M, 

whence the equation for our sought curve will be 

dgjf + gp) _ , P + 9{f + gp) 
z z ■ 

For expanding this equation, for the sake of brevity let us put p + g{f + gp) = S, 
and we will have: 

dgjf + gp) _dS__ SdZ 

z ~ z zz 

or 

dg{f^gp) = dS-^-^. 
Therefore because dS = dp + dg{f + gp) + gd{f + gp), our equation will be 

SdZ 



= dp + gd{f + gp) 



Z 

But on the other hand it is: 

dZ _ pdp+{f + gp)d{f + gp) 

Z l+pp+(f + gp)2 ' 

which should be multiplied by S = p + g{f + gp)- Then multiplying by the 
denominator 1 +pp + (/ + gp)^ we will have: 



= dp+{g- fp)d{f + gp) - gpdp{f + gp) + dp{f + gp) 

or 

= dp+{g- fp)d{f + gp) + fdpif + gp), 
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which we can then write in this form: 

= dp{l + ff + gg) + {g- fp){df + pdg). 



4. Although this equation is simple enough, it is not clear however how it 
could be reduced to differentials of the first degree. But I have observed that 
the problem can be dealt with by the following substitution, namely: v = ', 

whence it will he p = , and now by differentiating we deduce 
_ -{ff + 99)dv + (1 + vv){fdg - gdf) 

{f + gv? 

Further, it will be 

„ . v{ff + 99) 

9- JP= F-^ , 

f + 9v 

and next 

, ^ fdf + gdg + v{gdf-fdg) 
df + pdg = — . 

./ +5^^ 

With these substituted in, the equation becomes: 

= ^dv{ff + gg){l + f,f + gg) + v{ff + gg)ifdf + gdg) 
+ (1 + vv){fdg - gdf) + iff + gg){fdg -gdf). 



5. In order to turn this into a simpler equation let us set ff + gg = hh, 
and it will be fdf + gdg = hdh, and next let |- = fc, so that it would become 
fdg — gdf = ffdk, and thus our equation is contracted into this form: 

= -hhdv{l + hh) + h^vdh +{l + hh + vv)ffdk. 

Also since g = fk, it will be //(I + kk) = hh and then // = j^^, from which 
we will have: 

dk 

= -dv(l + hh) + vhdh +(l + hh + vv) — ; 

1 + kk 



next by putting v = sy^T+hh this equation is reduced to this form: 

, //. , , ^ dk{l + ss) 
= -dsV{l + hh) + ^^^. 

Therefore now the quantity s can be exhibited separately from the others, as 

ds dk 



1 + ss {1 + kk)^l + hh' 
which seems to be the simplest form to which it can be brought in general. 
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6. While we have attempted to determine the two variables y and x by the 
single variable z, since all three are involved in equal measure in the calculation, 
it appeared to me to treat this whole question so that all the formulae involve 
in equal measure the three coordinates x, y, z; this idea will be put to good use 
in the following investigations that I will take up. 

Supplement 

7. Let the differential equation of a given surface be: pdx + qdy + rdz = 0, 
where p, q, r are functions of the coordinates x, y, z: whence for this equation to 
be possible, this condition must be satisfied!! 

pdq — qdp qdr — rdq rdp — pdr ^ 
dz dx dy 

With this done, the following equation for drawing the shortest line on this 
surface will be obtained, which involves the three coordinates x, y, z is equal 
measure: 

ddx{qdz ~ rdy) + ddy{rdx — pdz) + ddz{pdy — qdx) — 0. 

Or if we put for the sake of brevity: 

dyddz — dzddy — /, 
dzddx — dxddz — g, 
dxddy — dyddx = h, 

it will be fp + gq + hr = 0, then indeed on the other hand fdx + gdy + hdz = 0. 
Next if an element of the shortest curve is put — ds, it will be ds^ = dx^ + dy^ + 
dz"^, then indeed also 

dds qddz — rddy rddx — pddz pddy ~ qddx 
ds qdz — rdy rdx — pdz pdy — qdx 

Application to a spherical surface 

8. Let the equation for the surface be xdx + ydy + zdz = 0, so that we have 
here p = x,q — y,r — z, and the first equation for the shortest path will be the 
following: 

ddx{ydz — zdy) + ddy{zdx — xdz) + ddz{xdy — ydx) — 0, 

whose complete integral is thus ax + /3y + 72 = 0, which is apparent from the 
nature of the matter. The question is therefore reduced to how this integral can 
be worked out. 

■^Translator: cf. p. 101 of D. J. Struik, Outline of a history of differential geometry: I, Isis 
19 (1933), no. 1, 92-120 
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9. Now with the carUcr equation fx + gy + hz = 0, if for this equation we 
put n = ^f-=^f it will hedn = d- ^f-'^f and then 

^ yax—xay^ yax—xdy 

^j-j zddx — xddz {zdx — xdz) {yddx — xddy) 
ydx — xdy {ydx — xdyY ' 

or by expanding 

X 

dH = - — ; —-7r((dyddz — dzddy)x + (dzddx — dxddz)y + (dxddy — dyddx)z) 

[ydx — xdy)'' 

and with the /, g, h that have been introduced it will be 

{fx + gy + hz) 



[ydx — xdyY 



Since on the other hand it is fx + gy + hz = 0, it will be dU. = and hence 11 is 
a constant quantity, which if we put = A, the differential equation of the first 
degree 11 = can be thus expressed: 

° yax—xdy ^ 

A{ydx — xdy) = zdx — xdz, 
which when divided by xx will be integrable; for it would become 

— = -+ B or Ay-Bx-z = 
x x 

or with the constants switched 

ax + (3y -\- ^z = 0. 

Since this equation is clearly everywhere drawn from the center of the sphere, 
great circles arise on the surface of the sphere; whence it follows that all great 
circles are all the shortest paths which can be drawn on the surface of the sphere. 

10. Since in these calculations everything is typically reduced to a single 

variable, if for effecting this we put dy = tdx and dx = udx, taking dx as 
constant, the first equation will be as follows: 

dt{') — pu) + du{pt — q) =Q. 

And the equation for the surface will he p + qt + ru = 0; whence, since it then 
becomes p = —qt — ru, the former equation takes this form: 

dt{r + qtu + ruu) — du{q + rtu + qtt) = 0. 

Next it will be 

f = dx^ {tdu — udt) , g = —dx^du, h = dx^dt, 
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then indeed ds^ = dx'^{l + tt + uu) and finally 

dds tdt + udu qdu — rdt pdu pdt 
ds 1 + tt + uu qu — rt r — pu pt — q 

11. And if desired to introduce as a fourth variable the angle by putting 
dx = td(p, dy = udtp, dz = vdip, the equation for the surface will be 

pt + qu + rv = 0. 

Next for the letters /, g, h we will have 

/ = df'^ {udv — vdu) , 
g = dif^ {vdt — tdv) , 
h = d(p^{tdu — udt), 

thus it will then he ft + gu + hv = 0. The equation for the shortest path will 
be: 

fp + gQ + hr= p{udv — vdu) + q{vdt — tdv) + r{tdu — udt) = 0, 
and it would finally become ds"^ = dip^ {tt + uu + vv) and hence 

dds tdt + udu + vdv qdv — rdu rdt — pdv pdu — qdt 
ds tt + uu + vv qv — ru rt — pv pu — qt 

12. Since no way presents itself to us for integrating the general equation 
given above for drawing the shortest path on a surface, even though many cases 
can be given in which integration of the equation for the curve succeeds, it will 
be worthwhile to have expanded several of these here as a conclusion. 

13. Let us begin with the case in which one of the quantities p, q, r vanishes. 
In particular, if r = 0, then the equation for the surface will be pdx + qdy = 0, 
in which case therefore the surface is a cylinder, whose base is determined by 
the equation pdx + qdy. And putting in the equation for ^ the given r = 0, 
it would become ^ = whose integral is Ids = Idz + la and hence, taking 
numbers, ds = adz or 

dx^ + dy^ + dz^ = aadz^ or dx^ + dy^ = dz^{aa — 1). 

It will then be 

z\/aa — 1 = / ^ dx"^ + rfy^. 



where / \/ dx^ + dy^ expresses an element of the base of the curve; whence it is 
apparent that the altitude 2; of a cylinder is always proportional to the arc of 
the base. 
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14. Let us consider the case in which p = x and q = y, where the equation 
for the surface will thus be: 

xdx + ydx + rdz = 0, 

in which all spherical bodies or rotated surfaces are contained. Then it will 
further be 

dds xddy — yddx 
ds xdy — ydx ' 

whose integral is 

Ids = l{xdy — ydx) + la, 
and hence by taking numbers 

ds dx^ + dy^ + dz^ 2 

— = xdy — ydx and so = [xdy — ydx) 

a aa 

or switching the constants 

(dx^ + dy'^ + dz'^)AA = {xdy - ydx)"^. 

For integrating this equation again let us put x = v cos (fi and y = vsin ip, and 
it will be 

dx^ + dy^ = dv^ + vvdip^; 

then indeed for the curve it will be vdv + rdz = 0, where r is some function of 
V, so that dz = - Next indeed it will be 

xdy — ydx = vvd(p; 

and with all these substitutions our equation will be: 

AA{dv'^ + vvdJ^ H ) = AAds^ = v'^dip'^, 

rr 

from which in turn we gather 

^ 2 _ AAjdv'^ + ^) _ AAdv''{rr + vv) 
— vvAA rrvv{vv — A A) 

and hence 

Adv I rr + vv 

dip 



I \ vv — AA 



15. For other cases our general equation above can be put to more use. 
And first indeed, because so far things have depended on a ratio between the 
quantities p, q, r, one of which will be assumed at our pleasure. Thus let r = — 1, 
so it becomes dz = pdx + qdy, and let p and q be functions of x and y with it 
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being (^) = (^)- Next let us put dy = irdx, and it will be dz = (p + TTq)dx. 
Then, taking the element dx as constant, so that it would be ddx = 0, it will be 

ddy = dndx and ddz = {dp + ndq + qd-K)dx. 

Now the three letters f,g,h can be expressed in the following way: 

/ = dx^{Trdp — pdn + TTTrdq), 
g = —dx'^{dp + -Kdq + qd'n), 
h = dndx^ . 

Indeed we have seen that the equation for the shortest path is 

Pf + 9q + hr = 0, 
which therefore takes this form: 

— d7r(l +PP + qq) + dp{np — q) + 'JTdq{np — q) = 

or 

d'iT{l +PP + qq) + {dp + 7rdq){q — np) = 0. 

16. Since the two formulae p + nq adn q — irp are the central components of 
this equation, it will be very helpful to work out a relation between them. To 
this end let us put ^"^^ = v, whence it is now tt = then indeed in turn 

v{pp + qq) 

g-Tip= . , 

p + vq 

and on the other hand it will be 

, , , pdp + qdq + v{qdp - pdq) 
dp + wdq = . 

p + qv 

Now if we put q = up, it will be 

u—v , , , du{l + vv) ~ dv{l + uu) 

TT = and then aw = . 

1 + uv (1 + uv)'' 

Next let us put pp + qq = tt, and since we let q — up, it will be 

t J J P J pdq-qdp 

pp = and a - - = du= , 

1 + uu q pp 

and hence 

, , , ttdu 

pdq — qdp = ppdu ■ 



1 + uu 

and when these values are substituted, because 

/'/,/, , , , td.t - ivlldu) : (1 + uu) 

'i~^P=~n~, \ dp + ndq= — r , 

p{l + vu) p{l + uv) 
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it will be 

^_ {du{l + vv) - dv{l + uu)){l + tt) vtt{t{l + uu)dt - vttdu) 
{1+uvY pp{l + uu){l + uvY 

or 

(1 + tt){du{l + vv) - dv{l + uu)) + vt{{l + uu)dt - vtdu) = 0, 
which can then be reduced to this form: 

du{{l + vv){l + tt) - vvtt) - dv{l + tt){l + uu) + vtdt{l + uu) = 

or this even more neat one: 

du 



1 + uu 



{1+VV + tt) - dv{l + tt) + vtdt = 0. 



Now let us put V = w^/T+ti, and it will be 

dv(l + tt) - vtdt 

dw = ^ 

{l+tt)i 

that is, it will be 

dv{l + tt) - vtdt = (1 + tt)idw; 

then indeed it will be 

l+tt + vv = {1+ tt){l + ww), 
and with these values substituted our equation itself thus becomes 

(1 + tt)(l + ww) - (1 + tt)idw = 0, 
1 + uu ^ ^ ^ ^ ' 

hence by separating we obtain 



du dw\/l + tt 



1 + uu 1 + WW 

and consequently 

dw du 



1 + ww {1 + uu)^/TTti' 

which equation can always be intergated whenever t is a function of u or when- 
ever pp + qq is a function of ^ or g a function of p. 

17. It will further turn out, so that g is a function of p, first if z and y are 
thus determined by x and another new variable w, so that it would be y = Ax 
and z = Bx, with A and B being some functions of u). Therefore since we have 
put dz = pdx + qdy, it will be 

Bdx + xdB = pdx + qAdx + qxdA, 
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where the terms involving the differential should be compared separately, whence 
it will he p = B — Aq; and comparing separately the terms containing the quan- 
tity X it will he q = ^ and hence p = ^'^'^^^'^^ And thus p and q are functions 
of uj and hence tt = pp + qq and u = ^ will be functions of the single quantity 

oj and ^/l + tt will be a function of u. Because of this, the equation found 
above for the shortest path admits integration. Moreover, in this case, namely 
in which y = Ax and z = Bx, a conical surface follows constructed upon any 
base. 

18. The equation given above would then be integrable when we put y = 

Ax + C and z = Bx + D; for then it will be 

dz = pdx + qdy = Bdx + xdB + dD 

and because dy = Adx + xdA + dC, it will also be 

dz = pdx + qdy = pdx + Aqdx + xqdA + qdC 

and hence, comparing the members containing the quantity x with each other, 
and then indeed those which are endowed with the differential dx, it will be 

B =p + Aq and dB = qdA, 

then 

dB , BdA - AdB 

q=— and p= . 

Besides these indeed it should be dD = qdC = or the functions A, B, C, D 

should thus be compared that dAdD = dBdC. If this occurs, p and q will again 
be function of the single variable uj and hence yl + tt will also be a function 
of u, in which case too it is possible to define a shortest path. This case seems 
indeed to complete all the planar surfaces which can be explained in the plane. 
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